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Abstract: An integrated silicon nitride resonator is proposed as an ultra¬ 
compact source of bright single-mode quadrature squeezed light at 850 nm. 
Optical properties of the device are investigated and tailored through 
numerical simulations, with particular attention paid to loss associated 
with interfacing the device. An asymmetric double layer stack waveguide 
geometry with inverse vertical tapers is proposed for efficient and robust 
fibre-chip coupling, yielding a simulated total loss of -0.75 dB/facet. We 
assess the feasibility of the device through a full quantum noise analysis 
and derive the output squeezing spectrum for intra-cavity pump self-phase 
modulation. Subject to standard material loss and detection efficiencies, we 
find that the device holds promises for generating substantial quantum noise 
squeezing over a bandwidth exceeding 1 GHz. In the low-propagation loss 
regime, approximately -7 dB squeezing is predicted for a pump power of 
only 50 mW. 
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1. Introduction 


Over the last decades an immense progress in control and manipulation of quantum states of 
light and matter and their mutual interaction has facilitated generation of squeezed states of 
light in a broad selection of physical systems: from the pioneering works on four-wave mixing 
in atomic vapors |Tj and optical fibres (2j over commonly used workhorse systems based on 
spontaneous parametric down-conversion in -materials to recent demonstrations of pon- 
deromotive squeezing through optomechanical interactions 00- 

In the emergent field of optical quantum technologies, squeezed quantum states of light have 
become an essential resource for many continuous variable (cv) quantum information and com¬ 
munication protocols @. And the feasibility of exploiting such quantum-correlated states for 
ultra-sensitive measurements has been demonstrated in numerous implementations of quantum- 
enhanced sensing, pioneered in the context of Mech-Zehnder interferometry 00, and by now 
applied at both micro- and macroscopic scales @00. However, the size and operational com¬ 
plexity of quantum light sources remains a limitation for practical and industrial applications of 
optical quantum technologies and integration with electronics. Efficient pulsed Kerr-squeezing 
in standard polarization maintaining fibres 111 TJ) requires fibre lengths on the order of tens 
of meters and even monolithic squeezed light sources 1 13-15) employ bulk nonlinear crystals 
that are orders of magnitude larger than lithographically defined electronic integrated circuits. 
Generation of twin beam squeezing has been demonstrated in whispering gallery mode disk 
resonators GD and recently in resonant on-chip structures G3- bringing the development of a 
fully integrated opto-electronic quantum sensing device an important step closer. But so far, an 
integrated source of quantum correlated cv states remains to be demonstrated. 

One major limitation for the integration of quantum light sources is optical loss. Free-space 
setups mainly suffer from Fresnel losses at optical interfaces, which can be mitigated by tai¬ 
lored anti-reflection coatings, reducing the effect to the 0.1% level per interface. However, in 
integrated waveguide circuits the light is also subject to material losses and scattering losses, 
originating from the roughness of the lithographically defined sidewalls and increasing quadrat- 
ically with the refractive index contrast An = n core — n c i ai jciing | fl8) . The former contribution is 
governed by intrinsic properties of the material’s band structure and cannot be compensated for. 
Bulk crystalline silicon has a bandgap of 1.1 eV leading to significant two-photon absorption 
(TPA) at wavelengths below roughly 2000 nm, and this is the reason why so much attention 
is currently drawn towards identification of CMOS-compatible alternatives to the otherwise 
so appealing silicon-on-insulator (SOI) platform for nonlinear optics. Scattering losses, on the 
other hand, can be compensated for by optimization of fabrication procedures and engineering 
of waveguide dimensions. In stoichiometric silicon nitride films deposited with low-pressure 
chemical vapor deposition (LPCVD), 200 nm thick channel waveguides with < —0.2 dB/cm 
propagation loss at 780nm have been demonstrated |18j|. Using temperature-cycling to over¬ 
come the tensile-stress limited film thickness of around 250 nm, propagation losses of -0.12 







dB/cm at 1540 nm have been measured in channel waveguides of more than 700 nm thick¬ 
ness G3- And by employing a large aspect-ratio waveguide geometry the dominating scattering 
loss contribution can be shifted from the sidewalls to the waveguide top and bottom surfaces, 
which for LPCVD Si 3 N 4 have roughnesses in the sub-nanometer regime, resulting in ultra-low 
waveguide propagation losses below -0.1 dB/m at 1580nm [20) . Unfortunately, this comes at 
the price of low field confinement, which is essential for efficient nonlinear processes. An al¬ 
ternative CMOS-compatible material is high-index doped silica (Hydex) which has recently 
been demonstrated to be a promising candidate for integrated nonlinear optics ©■ Due to its 
high linear refractive index (n = 1.5 - 1.9) and low propagation loss of ~ —0.06dB/cm a non¬ 
linear parameter of y ~ 233 W 'km 1 at 1550 nm is achievable in high-confinement channel 
waveguides [ 22j . 

Yet another severe loss-related challenge for high-index contrast waveguides is the inevitable 
mode conversion and effective index step associated with chip interfacing. The former leads to 
loss through coupling to radiation modes and the latter to Fresnel loss and low collection ef¬ 
ficiency because of the large numerical aperture of high-index contrast waveguides. The most 
common approach to overcome these complications is to employ inverse tapers (spot size con¬ 
verters) at the waveguide ends allowing near-adiabatic conversion between the tightly confined 
waveguide mode and that of a single mode tapered fibre or fibre pigtail. Using such coupling 
strategies approximately -1.5 dB/facet coupling loss has been reported for pigtailed high-index 
contrast Hydex waveguides [22) . 

In this paper, we study theoretically the feasibility of an integrated silicon nitride waveguide 
resonator as a source of quadrature squeezed light around 850 nm, conveniently centered in 
the Ti:sapphire tuning range. Classical optical parametric oscillations [23 ] and frequency comb 
generation [24 ) have been demonstrated in similar structures at telecom wavelengths, requiring 
careful engineering of dispersion properties in order to achieve a broad phase matching band¬ 
width. For the application discussed here, the phasematching condition is relaxed because the 
quantum correlated fields of interest are co-resonating rf-sidebands symmetrically distributed 
around the pump. This, however, introduces further complications at the state-interrogation 
stage, since the sidebands cannot easily be separated from the bright pump field and the phase- 
space rotations required for standard homodyne state tomography are not readily realizable. 
Integrated silicon nitride ring resonators have also been proven viable sources of photon pairs 
with controlable degree of correlation [25] . But so far, the cv quantum noise properties of 
the parametrically generated sideband fields have remained unexplored experimentally. Optical 
losses are detrimental to quantum correlated states in general but quadrature squeezed states 
are particularly vulnerable. The effect of optical loss on such states is not just a mere reduction 
in production rate but an over all degradation of quantum correlations and thereby it’s practical 
applicability. For a given amount of loss the state degradation is more pronounced for larger 
degrees of squeezing, rendering moderate squeezing levels a more practically tractable target. 

The primary objective of our analysis is to identify a feasible design for such an integrated 
device, subject to the constraints of conventional microfabrication techniques and commonly 
achievable propagation loss figures. We address this task by performing numerical simulations 
of optical key properties of the waveguide structure and using the results as inputs to a full 
quantum noise analysis of the nonlinear intra-cavity interaction. Special attention is paid to 
optimization of chip in and out coupling since this is in general a bottleneck for the integration 
of quantum light sources. 

2. Silicon nitride 

Silicon nitride is a CMOS compatible material widely used in microelectronics, both as gate 
dielectric in thin-film transistors and as mask, due to its high etching selectivity over silicon. 


In the field of integrated optics LPCVD deposited amorphous stoichiometric silicon nitride 
(Si 3 N 4 ) has only recently been introduced as a platform for nonlinear optics but has already 
been intensively investigated. Due to its linear refractive index of n ss 2 it is a suitable core 
material for high-index contrast waveguides in applications where tight field confinement and 
small waveguide bending radii are required. With a bandgap of approximately 5eV |26|, the 
absorption edge of silicon nitride is just below 300 nm, providing low material loss in the visible 
and infrared, and above 600 nm TPA does not represent a limitation even at high power levels, 
as opposed to silicon. For the particular application at hand, an equally important property of 
silicon nitride is the reasonably high third order Kerr nonlinearity 112 = 2.5 x 10 - 15 cm 2 /W 1271, 
which is essential for implementation of parametric processes, enabling manipulation of the 
quantum properties of light. For comparison, the nonlinearity of silicon nitride is roughly 10 
times the nonlinearity of silica, 2 times that of Hydex, and 0.1 times that of silicon. Thus, 
with sufficient optical isolation from the silicon substrate, e.g. using silicon dioxide as cladding 
material, silicon nitride is a promising candidate for integration of quantum optical circuits 
where minimization of losses and strong interactions are of uttermost importance. 


3. The system 

The specific integrated system that we are concerned with in this paper is a buried channel 
waveguide circuit consisting of a racetrack resonator (RTR) with radius of curvature R, laterally 
coupled to a straight bus waveguide (Fig. [TJ. A finite overlap of the evanescent fields of the 
two waveguide modes in the coupling region allows coupling to the RTR with an efficiency 
modelled by the field amplitude coupling rate y c ss k 2 /2t , where z = 2 « e ff(T c + nR)/c is the 
resonator round-trip time and K -2 is the intensity transmittivity of the coupling mirror of an 
equivalent one-sided bulk-optics cavity. The approximate expression for y c is valid for K~ <C 1. 
The coupling rate is controlled through the gap size and coupling length parameters g and L c , 
respectively. Intra-cavity losses, primarily due to scattering from the waveguide sidewalls, are 
represented by an intrinsic loss rate Yo ss /c ( 2 /2z = ac/ 2u e ff. where a is the waveguide per meter 
loss parameter, c is the vacuum speed of light, and « e ff the effective mode index. The total loss 
rate is given by 7 = 70 + 7c (HWHM) determining the loaded quality factor of the resonator. 
An important parameter for controlling the squeezed light generation efficiency of the system 
is the escape efficiency rjesc = 7c/ 7 , characterizing the collection efficiency of the intra-cavity 
field. The escape efficiency determines the coupling regime with rj esc = 1/2 corresponding to 
critical coupling and rj est < 1/2 (rjesc > 1 /2) to under (over) coupling. For the application of 
squeezed light generation, operation in the over coupled regime is required, but obviously there 
is a trade-off between efficient collection from the cavity (large rj esc ) and resonant enhancement 
of the nonlinear process (small 7 ). 

As depicted in Fig. [T] b), we opt for a rectangular low-aspect ratio cross sectional geometry of 
the waveguides. For the sake of minimizing fabricational complexity, we restrict the thicknesss 
t to a value of 250 nm, yielding the largest possible mode confinement while at the same time 
eliminating the need for stress-releasing temperature cycling in the deposition process. This in 
turn puts an upper bound on the width w, due to the requirement of single mode operation (cf. 
section | 6 T| ). 

4. Kerr nonlinearity 

The underlying mechanism for nonlinear optical effects in Kerr media is the third order tensor 

(3) 

susceptibility Xijki* supporting a great variety of four-photon mixing processes, collectively 
termed four-wave mixing (FWM). In general, FWM processes involve conversion of two pump 
photons fojp) into a pair of photons at signal and idler frequencies ( 0 J V , ft),), respectively, with 
the convention that co s > ft);. The interacting field frequencies are related by 2ft) p = ft) s + ft);. 
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Fig. 1. (a) Top view of the racetrack resonator geometry, (b) Cross sectional view at the 
bus-resonator coupling region. The thickness is t = 250 nm and the width w is chosen such 
that the waveguide only supports a single transversal mode. The gap size g is fixed, and 
limited from below by the resolution of the lithographic method chosen, and the coupling 
rate y c is determined by the length of the coupling region L c . 


required by energy conservation, and in order for the process to occur efficienctly the phase 
matching condition A/3 = 2/3 p - /3 S /3; « 0 must be fulfilled; here expressed in terms of the 

guided mode propagation constants /3,- = ;; e n ( OJ/)^o / where n e g- is the effective mode index. For 
this reason, the dominant interaction is the intrinsically phase matched self-phase modulation 
(SPM) process corresponding to the case of complete degeneracy of the partaking fields and 
resulting in an intensity dependent refractive index n(I) = «o +«2^- In this case the pump field 
experiences a self-induced phase shift (j ) n i = ImiilL/X = y n iLP p , where y n i = con 2 /cA e f{ is the 
nonlinear parameter, A e g the effective mode area, and P p the pump power. It is well-known that 
this type of FWM is capable of generating self-induced single-mode squeezing 128 29) and it 
is also the interaction of primary interest to this work. However, in the context of nonlinear 
optics in integrated resonators a more commonly discussed interaction is non-degenerate four- 
wave mixing (NDFWM) where an intense pump field interacts with two or a multiple of weaker 
non-degenerate signal and idler fields, each of the fields being resonant on separate longitudinal 
resonator modes. Parametric oscillators 1 23|30| and Kerr-frequency combs 1311 have previously 
been demonstrated using this interaction, and generation of squeezing, both below and above 
threshold, is discussed in |32) . 

A simple and intuitive understanding of pump squeezing by SPM is provided in (28) : Con¬ 
sider an input pump field a = (Xq + 8X\ n + 8Xf consisting of a coherent carrier amplitude 
and vacuum fluctuations in the sideband quadratures defined as SXq = 8ae lB + 8a' e‘ e with 
8X i = 8Xq _o and 8Xi = 5X e=n i 2 . Both the coherent amplitude ocq and the noise 8X'f drives 
the SPM process, resulting in a phase shift of the carrier proportional to ao and induced phase 
fluctuations at sideband frequencies given by <S@ n | = AKmLa.^SX'"/X. For small angles the cor¬ 
responding change in the phase quadrature fluctuations is A(<SX 2 ) ~ Oo<5</> n i = 2r8X'f, yielding 
the following input-output relations for the quadrature fluctuations under influence of SPM: 


8X° ut 


= 8X\ n 
SX 2 om = SXf + 2rSxl n 


(1) 

( 2 ) 


While the amplitude quadrature is unaffected by the interaction, an amount of the ampli¬ 
tude fluctuations proportional to the carrier intensity is admixed to the phase quadrature 
thereby correlating the two. The correlations are strongest at a quadrature phase of 0 m j n (r) = 
2 arctan(— 1/r) for which the variance is Var(5Xg Ut ) = 1 — 2r\J\ +r 2 + 2r 2 . Thus, for any r> 0 
the noise is, in principle, squeezed below the vacuum level at the optimal phase angle which in 




























the limit of weak interaction strength is lirrv >o 0 m in (f) - —n/4. However, the efficiency of the 
third order nonlinearity is very low and appreciable r-values are only achievable for very strong 
pump fields. In case of cw pumping this requirement can be addressed by resonant enhancement 
of the interaction which will be our focus in the following section. 


5. Quantum dynamics of resonantly enhanced SPM 


In this section we study the effects of SPM on the quantum dynamics and noise properties 
of a cavity field interacting with a Kerr nonlinearity. Quantum mechanically the interaction is 
governed by the quartic term 

fl ^ |2 2 

#SPM = ~Y a p (3) 

of the full FWM Hamiltonian (see Appendix [A]), where the interaction strength is given by 
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ficoc 2 y n i 

2 n efiL 


(4) 


Following the quantum Langevin approach we formulate linearized equations of motion for the 
intracavity quantum fluctuations, and the single mode squeezing spectrum for the output field 
is derived. 


5.1. Langevin equation of motion 


Considering a single pumped resonator mode a v interacting with a Kerr nonlinearity the general 
Langevin equation of motion is given by 133 341 


da p 

dt 


i 

h 


[flp,//] ya p + x /2y c a pjn e L 4“ 


(5) 


where 


H = h(O p a p a v + //spm- 


( 6 ) 


Operators fl p .; n and b v represent driving and input vacuum fields, respectively. Transforming all 
mode operators to a frame rotating at the driving laser frequency ft)| , by means of the substitu¬ 
tion fl p —> e~ l< ° Lt Ap, we get the rotating frame equation of motion: 


-rf = -(7- «'Ap)A p - i^A^A 2 + y/lf c A vM + ^/2^B pi 


(7) 


with detuning A p = 0 )l — ®p relative to the empty cavity mode frequency. Adapting the standard 
coupled-mode theory normalization convention, |A p | 2 represents the total mode photon number 
stored in the resonator and |A p ,j n | 2 represents the traveling mode photon flux. 


5.2. Linearized dynamics 

In order to solve the system dynamics we linearize the rotating frame quantum Langevin equa¬ 
tion Eq. 0, about the solution a p to the steady-state equation 

I 0£ p| 2 (T 2 + (A P —^|d£p| 2 ) 2 ) =2y c |ap )in | 2 . (8) 


We do so by making the substitution A p —> a p 4- 5a p and retaining only terms of first order in the 
quantum fluctuation operator 8a p . Defining the vector fluctuation operators (5a, = (Sa,. Sa- ) T 
the linearized system of differential equations can be stated as a matrix equation: 


t/(5ap 


dr 


[M - yl] <5a p + + ^/2]i ) 8b. 


(9) 







where the dynamics of the pumped mode is fully characterized by the system matrix 


M-yI = 


-r-K 2|e|-Ap) 

is* 


-i£ \ 

-Y+K 2|e| — A p ) J ’ 


( 10 ) 


and its eigenvalues 

A± = -y± \J\e\ 2 - (A p -2|e|) 2 . (11) 

For brevity we have introduced the pump parameter e = £ |ctp| 2 <?' 2 ^, where 0 is the phase of the 
pump field. From Eq. ( fT0| > we observe that SPM induces phase-dependent correlations in the 
fluctuations of the intra-cavity field (anti-diagonal elements), increasing with pump power, and 
a power-dependent nonlinear shift of the cold-cavity resonances (diagonal elements), responsi¬ 
ble for optical bistability of the system. 

Stability of the pumped mode requires that Re(A) < 0 for all eigenvalues. This is trivially 
fulfilled for A_, and for A+ the condition can be equivalently formulated as: 3|e| 2 — 4A p |e| + 
Ap + y 2 > 0. The real part of A + is plotted in Fig. |2 Solving the equation for |e| we find solu- 
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Fig. 2. Real part of the pumped mode eigenvalue A+ as a function of pump parameter and 
detuning. The A + = 0 contour (black dashed line) divides the parameter space into a stable 
(A+ < 0) and an unstable (A+ > 0) region. For detunings above the bistability condition 
(white dashed line) two stable regions exist and the system can be tuned between the two 
by changing the intracavity power. The intermediate unstable region is not accessible. By 
progressively detuning the drive field from the empty cavity resonance as the pump power 
is increased it is possible to maintain stability of the system and maximize the intra-cavity 
power (red dashed line). Furthermore, we observe that in terms of power the instability 
region is bounded from below by the condition that |e| > y indicating that the nonlinear 
scattering rate from the pumped mode should exceed the cavity dissipation rate in order for 
instability to occur. 

tions |e|± = j(2A p ± yj A 2 — 3y 2 ); for |e|_<|e|<|e|+ the system is unstable and outside it is 
always stable. Furthermore, for A p < y/3y no real positive solution for |e| exist and the stability 












criterion is always fulfilled; for A p > \/3y two solutions exist, marking the onset of optical bista¬ 
bility. We observe from Eq. (Ill that for the particular choice of detuning A p = t, |Op | 2 = |e| the 
eigenvalues collapse to a degenerate pair A± = — y. Thus, the system is always stable, indepen¬ 
dent of the pump power. In this case the SPM induced nonlinear shift of the cavity resonances 
is exactly compensated for by the pump laser detuning, restoring the usual proportionality be¬ 


tween input pump power and stored intra-cavity power, 
driven cavity. 




= 2 y c /y 2 • |ctp ; n |- for a resonantly 


5.3. Squeezing spectrum 

Since we are ultimately interested in deriving the squeezing spectrum for the system it is con¬ 
venient to transform Eq. (j9j) into frequency space. Fourier transformation of the differential 
equation system yields the following set of algebraic equations, 

- /£2<5a p (£2) = [M - yl] <5a p (£2) -{-^2y c <5a m (fl) + 5b(ft), (12) 

which can readily be solved for the intra-cavity fluctuations: 

<5a p (£2) = — [M - (y- (v^5a in (ft) + . (13) 

The intra-cavity field fluctuations can in turn be projected onto a set of output modes of the cav¬ 
ity using the input-output formalism. In this particular case the appropriate boundary condition 
is, 

Sa out (£2) = <5a; n (£2) - ^2^a p (£2). (14) 

Inserting Eq. ( fl3| l into Eq. ( fl4j i we find 

Sa out (£2) = [M— (y— i£2)I] -1 

x ([M + (Ay + i£2)I] 5 ain (£2) +2VWfc5b(£2)), (15) 

where Ay = y c — yo characterizes the coupling regime: Ay < 0 - under coupling. Ay = 0 - 
critical coupling, and Ay > 0 - over coupling. 

In general, an experimental interrogation of the cavity output field involves homodyne detec¬ 
tion of the field quadratures 

SX e (£2) = <Sfl 0Ut (£2)e“ ,e + 5a„ ut (-£2)e'' 0 , (16) 


for which the quantum noise properties are fully characterized by the normal ordered second 
order moment 


(:8Xe{Q)8X e (n!):) 


e 2l6 (8a out (^)8a out (^))+e 2,d (8al ut {Sl)8al ut {^')) 


+( Sa L t(-£2)5aout(£2')) 


<5(£2 + £V). 


(17) 


Following 133] the normal ordered spectrum of SXg (£2) is then found by integrating Eq. (17 1 
over £ 2 ', 

:Sq (£2) := f dQ!{: 8X e {Q.)8X g {Q !):), (18) 


and finally the squeezing spectrum of the output field quadrature is given by 


S 0 (£2) = 1 + :S 0 (£2): 


(19) 







which can be directly measured as the power spectral density of the homodyne photo current. 
The additional constant term accounts for the vacuum variance contribution from anti-normal 
ordered terms of the second order moment, which for the particular field quadrature definition 
in Eq. ( fl 6 | ) is equal to one. 

Using the SPM system matrix in Eq. ( p~0] > we can evaluate the contributing second order 
moments in Eq. (171. Combining terms results in the following expression for the squeezing 
spectrum of the pumped mode: 


nSPM 


(fl) = 1 + :Sf M (Q.): 


= 1 +G- 


2y|e| -2y[2|e| -A p ] cos2(0 + 0) 


where 


G = 


— [y 2 + fl 2 — Ap+4A p |e| — 3|e| 2 ] sin2(0 + 0) 


_ 4y c | £ | _ 

[Ajj + y 2 - Q, 2 - 4A p |e| + 31e| 2 ] 2 + 4y 2 fl 2 


( 20 ) 


( 21 ) 


Finite detection efficiency, primarily resulting from fibre-chip coupling (rj c ), imperfect mode 
overlap on the homodyne detector ( 77 mm ), and photodiode quantum efficiency (r/ qe ), causes 
admixing of vacuum noise and reduces the measurable squeezing according to: 


where // = t] c ■ J] n 


5 ™ eas (f 2 ) = (1 -17) +17 S g (£2) = 1 + 77 : S e (£ 2 ):, 
• r/qe is the total detection efficiency. 


( 22 ) 


6. Numerical simulations 

Up until now we have been considering the theoretical aspects of squeezed light generation in 
third order nonlinear materials, not accounting for actual material properties and system design. 
In this section we study the dependence of crucial optical propagation and nonlinear properties 
on geometrical system design parameters through numerical simulations. We fix the optical 
design wavelength to A = 850 nm, and the parameter space is constrained by aiming for a design 
with waveguide thickness t = 250nm and minimum lateral feature sizes >400nm. We do so 
in order to make the resulting design compatible with simple LPCVD single layer deposition 
and UV stepper lithography. All simulations have been carried out using the FieldDesigner and 
OptoDesigner software packages from PhoeniX Software. 


6.1. Waveguide mode properties 

A fundamental requirement for the waveguide architecture is that the transversal geometry 
should support only a single guided mode. This is important in order to optimize power con¬ 
finement in the waveguide core and to facilitate mode matching of the output field to a free 
space local oscillator field as required for homodyne interrogation. In Fig.[3]we plot the effec¬ 
tive mode indices for the first three guided modes as function of waveguide width w, derived 
from a full-vectorial mode analysis of the waveguide cross section. Starting from the square 
cross section and increasing the w: t aspect ratio the effective indices of the fundamental TE 
and TM modes separate with n™ > «™, indicating that the horizontally polarized TE mode 
is more confined to the silicon nitride core. This geometry-induced birefringence improves the 
polarization maintaining property of the waveguide and it is thus beneficial to maximize the 
aspect ratio. However, as the width is further increased the first higher-order TE mode starts to 
break through with the single-mode cutoff width being w co — 600 nm at the design wavelength. 






1.80 


1.75 


1.70 


u 1.65 

T3 

O 

g 

<D 

■3 1.60 

o 

w 

1.55 


1.50 


L45 —*—*—*— 1 —*—*—*— 1 —*—*—*— 1 —*—*—*— 1 —*—*—*— 1 —*—*—*— J — 

200 400 600 800 1000 1200 1400 

Waveguide width [nm] 



Fig. 3. Simulated (A = 850nm, t = 250 nm) effective mode indices for the first three guided 
modes as function of waveguide width. Insets: Simulated mode profiles corresponding to 
each of the five encircled points. Solid black lines outline the waveguide core cross section. 


An equally important property related to the waveguide cross section is the effective mode 
area A e ff introduced in section [4] Since the Hamiltonian coupling constant (Eq.|4| is inversely 
proportional to A e ff, minimization of this parameter will have a direct impact on the device effi¬ 
ciency. Different definitions of A e rr exist in the litterature, but here we adopt the form introduced 
by Rukhlenko et al. |35) given by 


^eff = «NL 


/ / Sy dxdy / S z dxdy , 
JJ-oo ' / J J NL 


(23) 


where anl is the cross sectional area of the nonlinear waveguide core and S, is the time- 
averaged z component of the Poynting vector, with the z-axis defined to be along the propa¬ 
gation direction of the bus waveguide. The effective mode area is thus simply defined as the 
core area scaled by the ratio of total mode power to power transmitted through the core. Simu¬ 
lated effective TEo mode areas as function of waveguide width are shown in Fig. |4]for a range 
of different thicknesses. 

Based on the simulations in Figs. 00 and the imposed fabricational constraints we con¬ 
clude that a cross sectional waveguide geometry with t = 250 nm and w = 500 nm would be 
a favourable choice meeting the requirements of being single mode, polarization maintaining, 
and having a small effective mode area. Further simulations will henceforth be based on these 
parameter values. 


6.2. RTR coupling efficiency 

We now turn to a discussion of the RTR to bus waveguide coupling efficiency and the im¬ 
portance of tailoring this parameter for the application of squeezed light generation. In many 
applications of micro resonators for nonlinear optics it is merely required that the circulating 
power be as large as possible, in order to maximize the nonlinear effects. In that case critical 
coupling is the targeted design. However, as briefly mentioned in section [3] it is as crucial for 
the application at hand that the cavity field predominantly decays through coupling to the bus 


















Fig. 4. Simulated (A = 850 nm) effective mode area as function of waveguide width for 
different waveguide thicknesses. Dashed lines show the nonlinear core area a^l which, 
according to the definition in ( |23| ), is the minimal effective mode area for a given cross 
section. 


waveguide and not via intracavity loss channels. The importance lies in the fact that any loss 
from the quantum correlated field is necessarily accompanied by admixing of vacuum fluctua¬ 
tions which will rapidly outbalance the induced quantum noise reduction and increase the noise 
floor back to the shot noise level. The figure of merit is the escape efficiency t] e sc = 7c/ 7 which 
should be as large as possible to ensure overcoupling and efficient routing of the cavity field, 
but at the same time a significant field enhancement factor of the cavity should be maintained. 

As mentioned, the intrinsic cavity loss rate 70 is primarily due to scattering losses, and for 
buried channel silicon nitride waveguides a conservative estimate for the propagation loss pa¬ 
rameter is a ~ -2dB/cm. A second contribution to the overall loss are bending losses due to 
sub-total internal reflection in waveguide bendings with small radii of curvature. To estimate 
the impact of this on the overall cavity loss, we have performed simulations of the bending 
loss suffered by the fundamental TE mode for different radii of curvature. The simulations are 
based on a combination of a 2D vectorial mode solver and coupled mode theory. Our findings 
are summarized in Table [T] and we observe that the effect of this loss channel can readily be 
reduced one order of magnitude below the linear propagation loss by increasing the resonator 
bending radius, while still maintaining a footprint of the structure well below 1 mm 2 . Finally, a 
third source of optical loss in RTRs are the transitions between straight and bending segments 
of the resonator. At these locations the optical field experiences abrupt changes in effective re¬ 
fractive index, acting as scattering centers within the waveguide. We will not pursue this loss 
mechanism further but just mention that the effect can be mitigated by bridging the transitions 
with adiabatically bending segments at the expense of an increased total resonator round trip 


length (36][37). 


Based on the above estimates the total round-trip loss is expected to be of the order of a 
few percent. To achieve strong overcoupling this has to be over-compensated for by the incou¬ 
pling efficiency, which in the case of an RTR geometry is controlled both through the gap size 
g and the length L c of the coupling region (see Fig. |T]». Figure [5] shows the simulated depen¬ 
dence of the coupling efficiency on g and L c . For large gaps the mode overlap is too small to 


























Table 1. Simulated (X = 850 nm, t = 250 nm, and w = 500 nm) bending loss for the TEo 
mode in resonators with different radii of curvature, R. 


R 

Simulated bending loss 

[j urn] 

(dB/360°] 

[dB/cm] 

25 

-1.62-10 3 

-0.103 

50 

-3.07-10- 3 

-0.0977 

75 

-4.35-10~ 3 

-0.0923 

100 

-5.59-10- 3 

-0.0890 


achieve efficient in and out coupling to the resonator whereas for very small gaps the energy 
oscillates back and forth between bus waveguide and RTR as function of the coupling length. 
Furthermore, we observe that for L c = 0 flm (ring resonator) the coupling efficiency drops off 
exponentially with gap size. Relevant simulation runs complying with the constraint on lateral 
feature sizes are plotted in the righthand panel. It appears that for g = 500nm the coupling ef- 
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Fig. 5. Simulated (X = 850 nm, t = 250 nm, w = 500 nm) efficiency of TEq coupling be¬ 
tween RTR and bus waveguide as function of coupling length L c . (left) Simulation runs 
with gap sizes in the range 100— 800 nm. (right) Selected simulation runs plotted on linear 
scale. Gray dashed lines indicate the total intra cavity loss for an RTR with R = 50 fim in 
case of -1 dB/cm (lower) and -2 dB/cm (upper) propagation loss. 


hciency can be conveniently tuned across and above the full range of expected intracavity loss 
values (gray dashed lines), rendering this gap size a favourable choice for fabrication of over¬ 
coupled resonators. Reducing the gap by another 100 nm leads to a strong perturbation of the 
intrinsic RTR properties and the associated increase in rj esc comes at a price of a much reduced 
finesse. Denoting the total intra-cavity power loss by «£? we can write approximate expressions 
for the RTR finesse and escape efficiency. 




2k 




l?e 


K? + 2zf’ 


(24) 


valid for K -3 , 2z? <C 1 • lable [ 2 ] summarizes the evaluation of (24 1 using simulation data from 
Fig. 0 We observe that for g = 500nm it is indeed possible to achieve escape efficiencies above 
80% and maintain & > 100, provided that waveguide propagation losses can be reduced to 
about -1 dB/cm. 



































Table 2. Evaluation of RTR (R = 50 ftm) escape efficiency and finesse in case of -1 dB/cm 
(left) and -2 dB/cm (right) propagation loss. Italicized values indicate that the approxima¬ 
tions used in ( |24| ) are only partially justified. 


L c 

g = 400 nm 

g = 500 nm 

Lc 

g = 400 nm 

g = 500 nm 

[Atm] 

f]e sc 

& 

^7esc 


[Atm] 

^7esc 

& 

^7esc 

& 

4 

0.88 

96 

0.73 

207 

4 

0.79 

86 

0.59 

167 

6 

0.91 

70 

0.79 

163 

6 

0.84 

65 

0.66 

137 

10 

0.94 

42 

0.84 

124 

10 

0.90 

40 

0.73 

108 

14 

0.96 

28 

0.88 

90 

14 

0.93 

27 

0.79 

81 


6.3. Robust chip coupling 

The notoriously large coupling loss associated with integrated waveguides imposes further se¬ 
vere limitations on the collection efficiency of the device. In this section we consider a waveg¬ 
uide design based on the TriPleX™ waveguide technology developed by LioniX and we inves¬ 
tigate numerically to what extend the design can improve on chip interfacing. Implementation 
of non-standard designs will inevitably increase fabrication costs and difficulties considerably 
but it is paramount that the coupling loss issue is addressed thoroughly in order for integrated 
quantum optical devices to be of any practical relevance. 

The design in question is sketched in Fig. [6] and consists of an asymmetric silicon nitride 
double layer stack with inverse vertical tapering of the upper layer. The lower film is defined 
by the same lithographic process as the upper layer, but unlike the latter it extends to the chip 
facet and acts as a weakly guiding structure for the incident light field serving two purposes: (i) 
it reduces interface Fresnel losses to ~4% (t7Fresnei = 96%), which can in principle be further 
reduced by anti-reflection coating, and (ii) the mode field diameter of the guided thin film mode 
can be easily tailored to match that of an incoming mode delivered from a tapered lensed fibre 
by controlling the film thickness. Light is then transferred to the main high-confinement waveg¬ 
uide by means of the inverse vertical taper. The efficiency of this process relies on optimization 
of the intermediate oxide layer thickness and the tapering angle £1 being sufficiently shallow 
for the mode conversion to take place adiabatically. 





Fig. 6. (a) Longitudinal cross section of the proposed silicon nitride double layer stack 
and inverse taper for robust chip coupling. Silicon nitride is indicated in red and the sur¬ 
rounding oxide cladding in blue. The performance of the design depends on optimization 
of thin film thickness t t f, intermediate oxide layer thickness t ox , and the tapering length 
L,. (b) Transverse cross section of the double layer stack. The lower thin film mirrors the 
lithographically defined pattern of the upper main waveguide, but with the essential excep¬ 
tion that it extends all the way to the chip facets, acting as a weakly guiding structure in 
connection with both in- and out-coupling. 


Based on mode analysis simulations of the fundamental TE mode of the lower stripe waveg- 






























uide as function of film thickness, we plot in Fig. [TJi the resulting mode waist size along with 
exponential fits. The resulting overlap between the asymmetric waveguide mode and a symmet¬ 
ric incoming Gaussian mode with a waist of 1 jlm, consistent with lensed tapered fibre coupling, 
is evaluated and plotted in Fig. |7 ]d (squares). For comparison, we also evaluate the closed form 
expression for the overlap integral of two symmetric Gaussian beams (waists wq, W 2 ) [38) 


V = 


4\v\w\ 
(wj+wj) 2 ’ 


(25) 


with w 1 = 1 jlm and the fitted exponential functions for w x and w y substituted for wj, respec¬ 
tively. Our calculations suggest that an overlap efficiency of ^overlap > 98% is achievable for a 
film thickness of 60 nm. For comparison, butt coupling a standard single mode fibre, with mode 
field diameter of 5 jlm, to the untapered waveguide would result in an unfeasible mode overlap 
of only 2%. 
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Fig. 7. Simulated (A = 850nm, w = 500nm) mode waist sizes along major (X) and minor 
(Y) axes of the lower stripe waveguide as function of film thickness (a), and corresponding 
overlap efficiency with an incoming symmetric Gaussian mode with a waist of 1 jlm (b). 

(c) Simulated delineation curves for approximate taper adiabaticity. Dashed horizontal lines 
correspond to linear constant-slope tapers of varying taper length. 

As can be seen from Fig. |7J) the mode overlap is a rather peaked function of film thickness, 
because of the exponentially increasing mode waist, meaning that the film thickness should be 
well controlled during fabrication. Using only a single silicon nitride layer the etching process 
defining the inverse taper must be terminated exactly at the targeted remaining film thickness, 
requiring a homogeneous and carefully calibrated etch rate. Furthermore, the etching process 
would result in considerable surface roughness of the thin film top surface causing increased 


































scattering losses. The proposed double layer stack design conveniently circumvents these fabri- 
cational issues: the thin film thickness is well-controlled through the LPCVD deposition process 
and the constraint on etch depth control can be relaxed because of the protective intermediate 
oxide layer. 

We now turn to the transfer of power from the thin lower stripe to the main waveguide. 
In order to put contraints on the waveguide taper design we adopt the theory developed by 
Love et al. f39]|40| for delineation of adiabaticity criteria for tapered optical fibres. Follow¬ 
ing this treatment, approximately adiabatic propagation of the fundamental waveguide mode 
along the taper is ensured if the local taper length-scale Zt is much larger than the coupling 
length Zb = 2 jt/(J3i — Pi) between the fundamental and second local modes everywhere along 
the taper. Here Pi = H e ff(®i)£o.i are the propagation constants for the involved modes, and for 
tapering angles -C 1 the local taper length-scale is zt ~ p /£2, where p is a local characteristic 
transverse length scale. For the linear constant-slope tapers cosidered here, zt is always equal 
to the position z along the taper. The condition z. t = z\, leads to the approximate delineation 
criterion 


Q.= 


P(P i ~Pi) 
27r 


(26) 


For angles much larger than Q. a significant amount of power is coupled to the second local 
mode and the taper becomes lossy. Figure |7J: shows delineation curves based on numerical 
mode analyses of the TEo and TMo modes along the taper. The adiabaticity criterion has been 
evaluated using both the local taper thickness fi oca i and the geometrical mean of the local cross 
sectional dimensions as characteristic length scale. Horizontal dashed lines indicate tapering 
angles corresponding to constant-slope tapers of length in the range 100 — 1000 pm. It is obvi¬ 
ous from the figure that any constant-slope taper will eventually enter the non-adiabatic region. 
However, by making the taper sufficiently long the associated loss can be reduced to a low 
level. In the lossy region approximately 10% of the power is coupled to the second mode per 
half coupling length, and as Zb increases dramatically as the taper thins down we estimate that 
the coupling loss can be reduced to < 10% for realistic taper lengths of 500 — 600 pm. To con¬ 
solidate this estimate, we have simulated the power transfer between lower thin film and main 
waveguide using 2D BPM simulations with the target parameter values derived in this section 
and using different intermediate oxide layer thicknesses. Without going further into details we 
mention that the simulation results suggest that a total transfer efficiency of the taper exceeding 
??taper = 90% is indeed realizable with an oxide thickness of t ox « lOOnm. 

In summary, the above analysis yields a promising estimate for the total per facet coupling 
efficiency for the asymmetric double layer stack design of 


Pc = PFresnel ' Poverlap ' Ptaper > 0.96 • 0.98 • 0.90 « 84%, (27) 

corresponding to -0.75 dB. Experimental confirmation of this value would be an important 
step towards feasible integrated high-index contrast structures for quadrature squeezed light 
generation. 

7. Source characterization 

The figure of merit characterizing the feasibility of the proposed integrated device is ultimately 
the maximum achievable quantum noise reduction. But other characteristica such as output 
state purity and bandwidth are also important for the practical applicability of the source and 
will as well be discussed in the following. 

In an earlier study of squeezing from cavity enhanced third-order nonlinearities by S. Rey- 
naud et al. |41j it was pointed out that perfect amplitude squeezing is, in principle, obtainable 
close to the turning points of the bistability curve (points on the black dashed contour in Fig. [2]). 



However, in the present study we will restrict our evaluation to the case where the pump field 
is adaptively kept on resonance with the hot cavity as the incident power in the bus waveg¬ 
uide is increased (dashed red line in Fig. [2ji. In this way the intra-cavity power is maximized 
and the steady state behaviour of the system is reduced to that of an empty cavity, deliberately 
eliminating optical bistability. 

The simulation results of section [6] provides a set of favorable target parameter values for 
the RTR design and the necessary inputs for estimating the resulting nonlinear interaction 
strength. Leaving the RTR power coupling efficiency as the only variable design parameter, 
we proceed to evaluate the theoretical quantum noise spectrum (20 p2| ) as function of exper¬ 
imentally accessible parameters external to the integrated device, yielding realistic values for 
key characteristics of the device. We assume a 250 x 500 nm waveguide cross section for which 
Fig. [d] yields an effective mode area of ,4 e fr = 0.1 8 pm 2 corresponding to a nonlinear param¬ 
eter of y n i = 10.2 W ! m . In the following, we further assume linear propagation losses to 
be -1 dB/cm and an RTR coupling gap of 500 nm. In order to accommodate the simultaneous 
requirements of high escape efficiency and appreciable intra-cavity field enhancement, we re¬ 
strict the coupling efficiency to values in the range 2-8%. Since the associated modification to 
the nonlinear interaction strength q. through the L c dependence of the RTR round-trip length, 
is only on the few-percent level, we will neglect this and simply assume a constant value given 
by the corresponding value for L c = 12 pm (k 2 = 5%), namely q = 117Hz. 

In order to estimate the measurable degree of squeezing the finite detection efficiency (22 1 of 
the setup must also be accounted for. In the previous section, a simulation based estimate for the 
per facet fibre-chip coupling efficiency was found to be T] c « 84%, and experimentally realistic 
values for the homodyne mode overlap and detector quantum efficiency are rj mm = T] qe = 98%. 
In total this amounts to a detection efficiency of T] = 80.6% which will be assumed in the 
following unless otherwise stated. 


7.1. Quantum noise tomography 

Using ( p0p2| i we calculate the quantum noise characteristics of the output field, accounting for 
imperfect collection and detection efficiencies. The corresponding relative noise power tomog¬ 
raphy is plotted in Fig.[8]as function of pump power for increasing escape efficiencies, assuming 
a particular sideband frequency of fl/27T = 30MHz well above any technical noise in a stan¬ 
dard experimental setting. Most significantly, a quantum noise reduction of about -4.5 dB is 
observable in the output field, confirming the feasibility of the designed device as a source of 
bright squeezing. Characteristic to the Kerr effect, both squeezing efficiency and phase depend 
on pump power and in conformity with the input-output relations (|2| the amplitude quadrature 
noise is unaffected by the interaction. 

As previously discussed, it is generally of interest to maximize the escape efficiency of a 
resonant source to preserve squeezing in the output field. However, for the system in question, 
design optimization is not as straight forward due to the self-driven squeezing interaction and 
power dependence of the Kerr nonlinearity: increasing rj esc inevitably reduces the cavity finesse 
and thereby the intra-cavity power driving the squeezing process. We observe from Fig. [8]that 
for a given |e| (bottom axes) larger squeezing is indeed achievable as the escape efficiency is 
increased. But the required pump power in the bus waveguide (top axes) increases accordingly, 
rapidly approaching experimentally intractable levels. For t] esc > 90% we find that reaching -3 
dB of squeezing already requires about 200 mW input pump power. 


7.2. State purity 

From the tomography in Fig. [8] we note that the anti-squeezing rises steeply with pump power 
compared to the squeezing level. To further characterize this property, we evaluate the state pu- 
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Fig. 8. (A = 850nm, a = —1 dB/cm, l] = 80.6%, CI/2tz = 30MHz) Relative noise power 
in the output field of a resonantly pumped RTR (A p = |e|) as function of pump parameter 
|e|. Panels (a)-(d) correspond to RTR escape efficiencies of 75%, 80%, 85%, and 90% 
respectively. Shot noise level is indicated by the dashed contour and the solid black contour 
corresponds to —3dB quantum noise reduction, (e) Cross sectional views of panels (a)-(d) 
corresponding to an input pump power of 200 mW. 











































rity, which for Gaussian states is defined as fx = 1 / ^/Dct [cr] [42 , a being the covariance matrix 
for the state. Subject to optical loss, the purity of a squeezed state will drop quickly due to ther- 
malization with the vacuum, rendering it practically impossible to produce highly squeezed 
minimum uncertainty states saturating the Heisenberg inequality, Var(5X)) • Varl'dA'?) = 1. Due 
to the comparatively large intra-cavity and coupling loss associated with the RTR source, con¬ 
siderable excess anti-squeezing and correspondingly low state purity is anticipated even for 
modest squeezing levels. Using the simulated data in Fig.[8]we evaluate the state purity as func¬ 
tion of squeezing, plotted in Fig. [9} both for the estimated experimental detection efficiency and 
in the idealized case of unity efficiency. As expected, the purity increases with r] esc but at the 



Relative noise power [dB] 



Fig. 9. State purity as function of noise reduction for different RTR escape efficiencies. 
Dashed segments indicate that generation of the corresponding noise reduction requires 
a pump power in the bus waveguide exceeding 200 mW. The two panels correspond to 
detection efficiencies of rj = 80.6% (left) and rj = 100% (right). 


expense of large pump powers. To circumvent this, a lower intra-cavity loss rate is required. 
Furthermore, Fig.[9](right) shows that despite the efforts to minimize fibre-chip coupling losses 
the device collection efficiency remains a limiting issue. It is clear that the inherent loss chan¬ 
nels of the source pose a severe limitation to the performance which can only be overcome by 
improved fabricational quality of the device. For comparison we plot in Fig. 10 the performance 
of the device with waveguide propagation loss and collection efficiency set to -0.2dB/cm and 
95%, respectively. This allows for an increase in escape efficiency and a dramatic reduction of 
the required pump power, enabling -6 dB of squeezing with a purity of /x « 0.8 using less than 
50 mW pump power. 



Fig. 10. State purity as function of noise reduction with device loss properties of a = -0.2 
dB/cm and t] c = 95%. Dashed segments indicate that generation of the corresponding noise 
reduction requires a pump power in the bus waveguide exceeding 50 mW. 































7.3. Squeezing bandwidth 


Finally, we investigate the noise spectrum of the RTR output field. Commonly, cw parametric 
down-conversion sources of quadrature squeezed light are limited to squeezing bandwidths 
of a few tens of MHz, set by the OPO cavity bandwidth. One exception being monolithic 
cavity systems ( |43]|44) ) where the internal loss is sufficiently low to reach bandwidths up 
to 170 MHz. Weak squeezing with a bandwidth of more than 2 GHz at 1550 nm has also been 
demonsUated in such a system solely resonant for the pump field ]45) , and in a later work about 
2 dB of squeezing was observed with a bandwidth exceeding 1 GHz 1461. 

Squeezed light is a widely used resource in quantum communication applications, e.g. for 
establishing the essential two-partite entanglement for quantum key disttibution |[6 [47 1. In this 
particular case the key rate is directly linked to the squeezing bandwidth and demonstration 
of efficient GHz-bandwidth sources is thus an important prerequisite for high-speed protocols. 
Large-bandwidth squeezed states are also essential for extending the applicability of cv quan¬ 
tum sensing techniques. As an example, state of the art optomechanical photonic crystal devices 
have mechanical resonance frequencies in the GHz range [ )48| , rendering conventional squeezed 
light sources insufficient for enhancing the transduction sensitivity. The integrated source de¬ 
sign discussed in this paper has potential to fulfill those requirements. As shown in the preced¬ 
ing sections, efficient squeezing is indeed possible and due to the intrinsic phase matching of 
the SPM process the spectral extend of the squeezing is limited by the cavity bandwidth alone. 



Sideband frequency [fi/yj 


Fig. 11. (X = 850nm, a = — ldB/cm, |e|/y = 0.5) Spectrum of the squeezed and anti¬ 
squeezed quadratures assuming 80.6% (thick), 91.2% (thin), and unity (dashed) detection 
efficiency. Panels (a)-(d) correspond to RTR escape efficiencies of 75%, 80%, 85%, and 
90% respectively. 


Figure [TT] shows the spectra of the squeezed and anti-squeezed quadratures corresponding 
to |e|/y = 0.5 and for the same escape efficiencies as in Figs. [8] and [9] With a total detection 
efficiency of 80.6% (?] c = 84%) the output field shows more than -2 dB quantum noise reduction 
over a bandwidth larger than 1 GHz across all four rj eS c-values. If the detection efficiency is 
increased to 91.2% (rj f = 95%) we expect more than -3 dB squeezing over a bandwidth of 
1.5GHzfor?7esc >85%. 
























































8. Conclusion 


While on-chip photonics is already a promising technology gradually evolving into a versatile, 
modular, and reconfigurable architecture, development of an integrated platform for cv quantum 
optics has only received little attention. That despite the vast toolbox of Gaussian cv quantum 
sensing, computation, and communication protocols, already successfully implemented and 
demonstrated in free-space optical systems, but unsuitable for marrying current technology 
without a scalable architecture and cost-efficient fabrication line. Squeezed states are one of the 
most commonly deployed resources for quantum-enhanced sensing applications and entangled- 
state generation for communication protocols, and as such development and demonstration of 
an efficient integrated source of quadrature squeezed states is an important step towards a cv 
state-based quantum-enhanced optical technology. 

We have presented a theoretical feasibility study of a resonant integrated source of bright 
single-mode squeezed light, exploiting the third-order optical nonlinearity of silicon nitride. 
Through numerical simulations we have studied the dependence of optical propagation and 
nonlinear properties on waveguide cross section geometry as well as resonator coupling and 
chip interfacing. The issue of optical loss in integrated optics and its detrimental influence 
on optical quantum states has been addressed. To this end we have advocated an asymmetric 
double layer stack waveguide geometry with inverse vertical tapers as a means for mitigating 
the inherently large coupling losses associated with interfacing high-confinement waveguides. 
On basis of the simulations we have proposed a concrete design optimized for generation and 
collection of squeezed light, and the feasibility of the design was evaluated from the derived 
squeezing spectrum of the output field. The theoretical predictions are encouraging, suggesting 
that the source is capable of producing about -4.5 dB of squeezing with 200 mW pump power, 
mainly limited by the chip-fibre out-coupling efficiency and under the conservative assumption 
of -1 dB/cm propagation loss. Furthermore, a quantum noise reduction exceeding -2dB is ex¬ 
pected over a bandwidth larger than 1 GHz. If the waveguide propagation loss can be reduced to 
-0.2 dB/cm and the chip out-coupling efficiency increased to 95%, then squeezing efficiencies 
of about -7 dB are attainable using only 50 mW pump power. 

A. Derivation of the Kerr interaction Hamiltonian 

In this appendix we briefly review the derivation of the quantum mechanical Hamiltonian de¬ 
scribing Kerr interaction in a dielectric waveguide. We will limit our analysis to involve only a 
single spatial mode of the waveguide and linearly polarized monochromatic fields. 

Classically, the field energy in a dielectric medium is given by j f dV {$(r,t) ■ @(r,t)), 

where the electric displacement is given by @(r,t) = £o<^(r,f) + &(r,t) and the brackets de¬ 
note cycle average. In case of a nonlinear dielectric, the anharmonic response is accounted 
for by Taylor expansion of the polarization 3?(r,t) in orders of the susceptibility, and for 
a Kerr medium the lowest non-vanishing contribution is of third order, characterized by the 
sucecptibility tensor of rank 4. In general, calculation of the corresponding interaction en¬ 
ergy involves the contraction of the susceptibility tensor with four individual fields - two pump 
fields ( cipi , a p 2 ) plus upper (a + ) and lower (a ) sideband fields. However, in the case of linearly 
co-polarized fields the expression simplifies significantly: 

M = dV(S{r,tf), (28) 

where (r,f) is still assumed to be a superposition of four independent waves. 

Since we are concerned only with a single spatial mode of the waveguide, an appropriate 



plane-wave expansion of the quantized electric field inside the medium is given by 


E(r,t) = iS\(ate^-^+ale-^-wA 
i 


where i = {pi ,p 2 , +, -} (29) 


with normalization constant So = yJh(o/2£on 2 V , where the quantization volume is V = L-A e g- 
given by the product of the effective area of the guidede mode and the waveguide length. We 
have assumed that the frequencies of the interacting modes are sufficiently close that dispersion 
effects and the difference in photon energies can be neglected. Replacing S(r,t) in (28 1 by the 
quantized field operator results in a large number of quadrilinear products of the four involved 
mode operators, of which all the energy non-conserving terms are readily removed by the cycle 
average, enforcing the rotating wave approximation. The remaining terms are all permutations 


of the form a/Oyatat (and Hermitian conjugates) for which OJ, 


(Oj = a>k + ft)/ and the phase 


factor associated with these terms is e + '* k,tk ( k,) ' r = e +,Ak r = e +,Ak '\ where in the last step 

we have used the assumption of a single mode waveguide. Carrying out the volume integral 
yields 


[dVe iAkz =[ [ dxdy [ dze iAkz = A eff L<i>(AkL) 
J J J Apfe J 0 


(30) 


(or the complex conjugate) where we have defined the phase matching function <-l>(AkL) = 
e' AkL / 2 sinc(AkL/2). The effect of wavevector mismatch on energy transfer efficiency of the 
processes is entirely given by |<t>(ACL)| 2 = sinc 2 (AkL/2). The factor e ,AAi / 2 only constitutes 
an immaterial global phase and will be neglected. Using ( |30| ) and the expression for Sq we can 
express the pre-factor of the quantized Hamiltonian, characterizing the interaction strength E: 


2 = \^ (/s|W) (31) 

= (32) 

= /i<i;<f>(ACL) where B, =- ,^ iL . (33) 

In^L 


Here we have used the definition of the nonlinear refractive index «2 = 3/y 3 */4fc'o// 2 c and the 
nonlinear parameter = njOj/cA^f. 

We now return to the sum of operator products, but we restrict our treatment to the case of 
degenerate pump fields (a p \ = a p 2 = a p ). Furthermore, we choose to write the Hamiltonian in 
normal ordered form. Following the treatment by Jack et al. |32) the Hamiltonian is then given 
by 

Hi = 2^ : CL‘( a ‘+ a l)) 4: where i = { P)+j _} (34) 

and we are implicitly assuming that energy non-conserving terms are neglected. The overall 
factor of 2 accounts for anti-normal ordered terms in the expansion and the 4! in the denomina¬ 
tor is for normalization with respect to the total number of permutations. The remaining energy 
conserving terms and their interpretation are listed in Table |3]below. The interactions described 
by the three associated processes are sketched in Fig. |T2]in case of resonant interaction between 
a strong pump field and two weaker signal and idler fields, all resonant on separate longitudinal 
cavity modes. 

From the definition of <-l>(AkL) it is clear that only the FWM proces requires active phase 
matching whereas both SPM and CPM are intrinsically phase matched. The relative strength 










Fig. 12. Pictorial illustration of the x ^-mediated processes SPM, CPM, and NDFWM, 
in the case where a strongly pumped resonator mode at frequency (Up interacts with two 
weaker signal and idler sideband modes at frequency co s and fit;, respectively. 


Table 3. Energy conserving operator products in the expansion of the normal ordered 
interaction Hamiltonian ([34} 

Operator terms 

Associated physical process 

+2 9 f 2 2 t 2 2 

a ^ a + a+ a_ a_ 

2 t t 
a^a + a_ 

a^a v a^a~ a+a+ala- 

Self-phase modulation (SPM) 

Non-degenerate four-wave mixing (NDFWM) 
Cross-phase modulation (CPM) 


of the three processes can be evaluated by using the multinomial distribution for expanding the 
Hamiltonian. Calculating the multiplicities for the three processes we get weight factors of 1/2, 
1, and 2 for SPM, FWM, amd CPM, respectively. Finally, we can write the total Kerr interaction 
Hamiltonian (assuming perfect phase mathcing): 


/i^<f>(AkL)(apfl i j_at +a^ 2 a + a~) 

(35) 

+2h£ (a^a ? a\_a+ a^a^a^a- +a^a + ata_) 

(36) 

. +2 9 

+ —(ctp a v + a + a + +a'_ a_). 

(37) 
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